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It was proposed that the id^2_y2 density-wave state (DDW) may be responsible for the pseudogap behavior 
in the underdoped cuprates. Here we show that the admixture of a small dxy component to the DDW state 
breaks the symmetry between the counter-propagating orbital currents of the DDW state and, thus, violates the 
macroscopic time-reversal symmetry. This symmetry breaking results in a non-zero polar Kerr effect, which has 
recently been observed in the pseudogap phase. 
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Introduction. The nature of the pseudogap in the cuprate 
high-temperature superconductors is a much-debated and still 
unresolved issue |1]. Many of the anomalous properties as- 
sociated with the pseudogap can be best described as result- 
ing from the order parameters |2], such as those of various 
density-wave states, competing with the d-wave superconduc- 
tivity. If a density-wave state were to have an energy gap of the 
d-wave symmetry on the Fermi surface, in conformity with 
the pseudogap, the natural candidate would be the idx2_y2 
density wave (DDW) state. Indeed, much of the phenomenol- 
ogyof the cuprates in the underdoped regime can be unified 
||3l|4j,pj by making a single assumption that the ordered DDW 
state is responsible for the pseudogap. However, the situation 
is still controversial, because there is no convincing direct ex- 
perimental evidence for a phase transition to another state in 
the underdoped regime, besides the antiferromagnetic order 
close to half-filling. 

The recent observation of quantum oscillations ^ in un- 
derdoped YBCO in high magnetic fields has renewed interest 
in the ordered DDW state [7]. The special coherence factors 
of the DDW state reduce the hole pockets to appear as Fermi 
arcs in angle resolved photo emission experiments An 
important recent development in the context of the pseudo- 
gap is the observation of a non-zero polar KeiT effect (PKE) 
in the underdoped YBCO [9], which demonstrates the macro- 
scopic time-reversal-symmetry breaking. The PKE appears 
roughly at the same temperature where the pseudogap devel- 
ops (isd, and, near the optimal doping, appears at a temperature 
below the superconducting transition temperature. This obser- 
vation suggests that the time-reversal symmetry breaking and 
the pseudogap in the cuprates have the same physical origin, 
which is unrelated to the superconductivity. Similar conclu- 
sion was reached earlier, using muon spin relaxation experi- 
ments, in Ref. jioll . In this Letter, we examine the PKE of the 
ordered DDW state, and show that it can give rise to a non- 
zero signal only if there is an admixture of a d^y component 
to the total order parameter amplitude. There can be different 
mechanisms for spontaneous generation of the d^y potential 
in the mean-field Hamiltonian due to microscopic interactions 
between electrons or due to a structural transition with 
the appropriate symmetry. It is important to stress, however. 



that the amplitude of the dxy admixture, as deduced from the 
PKE signal [jgt], is tiny relative to the main idx2_y2 compo- 
nent. Thus, it may be difficult to uncover it unless the ex- 
periment directly probes broken time-reversal symmetry. Su- 
perconducting order parameters of similar form, explaining 
experimental observation of broken time-reversal symmetry, 
have been discussed in lfl2ll . 

idxi^y2 and {d + id) density-wave states. The commen- 
surate DDW state fll'l is described by the following order pa- 
rameter, which is a particle-hole singlet in spin space, 

/ t \ 

\fk+Q,a^k,f}j o^iWkScf), Wk = — (cosfc^ -coaky). 

(1) 

Here and c are the electron creation and annihilation oper- 
ators on the square lattice of the copper atoms, k = (kx, ky) 
is the two-dimensional momentum, Q = (tt, tt) is the wave 
vector of the density wave, and a and (3 are the spin indices. 
Since the order parameter involves Sap, we do all of our cal- 
culations ignoring the spins and just multiply the final results 
by two. We also set the lattice constant a, the electron charge 
e, and the Planck constant h to unity (e, h,a = 1) at the in- 
termediate stages of the calculations and restore the full units 
only in the final results. In Eq. iWk is the DDW order pa- 
rameter with the idx2^y2 symmetry in the momentum space. 
For Q = (tt, tt), it is purely imaginary flT] and gives rise to 
spontaneous cuiTents along the bonds of the square lattice, as 
shown in the left panel of Fig. [T] Although the presence of 
the spontaneous currents in the DDW state violates the micro- 
scopic time-reversal symmetry, there is no macroscopic vio- 
lation. This happens because the DDW state preserves the 
combined symmetry of time reversal and translation by one 
lattice period. As a result, the staggered magnetic flux pro- 
duced by the cuiTents averages to zero on the macroscopic 
scale. So, one has to go further in order to explain the macro- 
scopic time-reversal symmetry breaking revealed by the PKE 
measurements |@] . 

We propose to consider a generalization to the d^y + 
idx2^y2 density wave, where the order parameter is a com- 
bination of two density waves with different angular patterns: 

(4+Q,a^fc,/3^ OC [iWk - Afe] Sap, (2) 
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FIG. 1: Left: The real-space picture of the id^2_y2 density- wave 
state. Spontaneous currents are shown by the arrows. Clockwise 
and counterclockwise circulations are indicated by the — and the + 
signs, respectively. Right: The same for the dxy + idx2_y2 density- 
wave state. Now the neighboring plaquettes have different ampli- 
tudes of the next-nearest-neighbor tunneling, indicated by the solid 
and dashed diagonal lines. 



Afc Aq sin kx sin fc„ 



(3) 



The real-space structure of this order parameter 01 111 is shown 
in the right panel of Fig. [T] The solid and dashed diagonal 
lines represent the staggered modulation of the next-nearest- 
neighbor electron tunneling due to the dxy component of the 
density wave. This staggered modulation breaks the symme- 
try between the plaquettes with positive and negative circula- 
tions, so the macroscopic time-reversal symmetry is violated. 

The Hartree-Fock Hamiltonian describing the mean-field 
{d + id)-DW is 



H = 



fj. iWk - Afc 
- Afc £fe+Q — M 
2t(cos kx + cos ky) + At' cos kx cos ky, 



£k - 



(4) 
(5) 



where is the band dispersion of the electrons, and fi is 
the chemical potential. The Hamiltonian (|4| operates on the 
two-component spinor 5'^ = (cfc,Cfc+Q) and can be ex- 
panded over the Pauli matrices r and the unity matrix /, 



H 



wo{k)I + w{k) 



£fc+gfc+Q 



The spectrum of the Hamiltonian consists of two branches 
with the eigenenergies E±{k) — wo{k) ± w{k), where 
w{k) — \w{k)\. For a generic set of parameters, corre- 
sponding to a non-zero hole doping in the underdoped regime, 
the reconstructed Fermi surface consists of two hole pock- 
ets near the {tt/2, ±7r/2) points and one electron pocket near 
the (tt, 0) point in the reduced Brillouin zone, as shown in 
Fig.|2] The Hamiltonian (|4]l with t' = and fi = was stud- 
ied in Ref. where it was found that this system exhibits 
the intrinsic quantum Hall effect without an external magnetic 
field. This confirms that the time-reversal symmetry is spon- 
taneously broken in the {d + id) density-wave state. 

Polar Kerr effect in the d + id density-wave state. The po- 
lar Kerr angle 6k can be expressed in terms of the imaginary 
part of the ac Hall conductivity cr"j,(rj) at a frequency il Ill4l] . 
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FIG. 2: The Fermi surface of the dxy + idx2_y2 density-wave state, 
consisting of the electron and hole pockets. The hole pockets are 
located near {n/2, ±7r/2), and the electron pocket near (tt, 0) in the 
reduced Brillouin zone. The band-structure parameters are given at 
the end of the paper, and Wo = 20 meV is taken for this plot. 



where axy is the Hall conductivity per 2D layer, which is re- 
lated to the 3D conductivity by the interlayer distance d, and 
n is the refractive index of the system. The Kerr ang le was 
calculated for a px + ipy chiral superconductor BlSll . moti- 
vated by a non-zero PKE in Sr2Rti04 [16]. Eq. ^ is valid 
when the imaginary part of n is small relative to the real part. 
Otherwise, there is also an additional contribution ^4] from 
the real part a'xyi^^)- For simplicity, we focus on (T"j,(f7), but 



a'xy{il.) can be also easily obtained from our calculations. 

Introducing the Lagrangian density, C = ioJnl — H{k), 
where iv^ is a fermionic Matsubara frequency, we obtain the 
electron Green's function as a function of fc = (ia;„, kx,ky) 



Go(fc) 
9ik) = 



woik)]' 



wo{k)] i + w{k) ■ T 

9{^n-,k) 
-\w{k)\\ 



(7) 
(8) 



In order to calculate the Hall conductivity, we consider the 
electromagnetic potential A = {Ao,Ax, Ay), where Aq is the 
scalar potential, and A = {Ax, Ay) is the vector potential. It 
is introduced into the Lagrangian density £ by the standard 
substitution k ^ k A and iujn iujn — Aq, where c is 
absorbed into A. Expanding £ to the first order in A, we find 
the following addition to the Lagrangian density 



r(g, k) = -Ao(q) i + A{q) ■ Vfc [w{k) 



(9) 



Here the electromagnetic potential A (q) is assumed to be a 
function of the Fourier variables q = {inm,qx,qy), where 
ri„i is a bosonic Matsubara frequency. The Lagrangian ^ 
operates between the fermionic spinors 'i>{k) and ^^(fc + q). 
Calculating the one-loop Feynman diagrams, we find the fol- 
lowing effective action to the second order in F, 

S = ^TTY,^{q,k)Go{k)T{^q,k + q)Go{k + q). (10) 

k,q 

Substituting Eqs. and Q into Eq. ( fTOl i. we write the 

effective action in the form 



k,q 



Tr Ci 

Co ' 



(11) 
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where Ci is the product of the numerators from Eq. (O and 
Eq. ^, and C2 = g{k + q) g{k). 

Calculating the trace over the f matrices in the numerator 
of Eq. (fTTT i. we look for the Chern-Simons-like terms, which 
are responsible for the Hall effect. Indeed, we find two such 
contributions to the effective action, 5*^ = 51 + 5*2, 

5*1 = i^<Ta;yAo{cf)[qxAy{~q) - QyAxi-q)], (12) 
^2 = X]^a.,[A,((?)A,(-9-)-A,(q-)A,(-9)].(13) 

By taking the trace in Eq. ( ITTI ). we find the following expres- 
sion for the coefficient axy in Eqs. ( fT2] i and ( fTST l 



i?(fc) 



3^ 5(u;„, fe) g{uJn + ilm, fc) ' 



B{k) = -2w 



dw dw 



(14) 



(15) 



where the vector function w{k) is taken from the Hamilto- 
nian. Indeed, the numerator in Eq. (fTTT i is a product of four 
terms, each containing the f matrices. The Chern-Simons ac- 
tion (fT2] i and (T3[ is obtained when these terms contribute dif- 
ferent f matrices and the unity matrix to the trace. Using the 
formula Tr[fi7Vn,f„] = 2ie;m„, where eimn is the antisym- 
metric tensor, we arrive at Eq. ( fTSl l. We also expanded the 
numerator of Eq. (fTTT i to the first order in q in order to obtain 
Eq. (fTSl i and then took the limit q 0, which is appropriate 
for the electromagnetic response of non-relativistic electrons. 
Substituting the expression for w{k) from the Hamiltonian 
into the general formula (fTSl l. we find B{k) for the d + id 
density wave, B{k) — 4iAoWo(sin^ ky + cos^ ky sin^ kx). 
Notice that B{k) is proportional to the product of Ag and Wq, 
so the intrinsic Hall conductivity is non-zero only when both 
components of the d + id density wave are present. The func- 
tion B{k) in Eq s. (ITsT i and above is the Berry curvature for 
the Hamiltonian lll7ll . 

The coefficient axy in Eqs. (fT2] | and (fT3] i is the Hall conduc- 
tivity. Indeed, taking the variation jx{q) — SSh /SAx{—q), 
we find the Hall relation at a finite frequency 



{n„,)Ey{n 



(16) 



where E = —iqAo — is the electric field. Now we 

perform summation over the fermionic Matsubara frequency 
u!n in Eq. ( fT4l l and analytic continuation to the real bosonic 
frequency iftm ^ + i6, where S is infinitesimal. The result 
is 



axy{^) 



Let us now explore different limiting cases for Eq. ( fTTI l. 
First we consider the dc limit il = at T = in the case when 
there are no pockets, i.e. when i?+(fc) > and E^{k) < 
for all k. This situation takes place, e.g., when t' = and 
/i = 0. In this case, using Eq. (fTsT i, we can write Eq. (fTTI l as 



RBZ 



d^fc B{k) {NF[E+{k)] ~ NpiE^jk)]} 
(27r)2 w{k) [n + i5- 2w{k)] [fl + iS + 2w{k)] 

(17) 

where Np{E) is the Fermi occupation function for a temper- 
ature T, and 2w{k) — E+{k) — E^{k) is the energy gap 
between the bands. The integral over k is taken over the RBZ. 
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where n(fc) = w{k)/w{k) is a unit vector. The integral 
in Eq. ( fTSl l is an integer topological invariant, which repre- 
sents the degree of mapping from the vector k to the unit 
vector n. Using Eqs. (fTSl l and ( fTSl l. we indeed find the in- 
teger quantum Hall effect per one spin component C7xy{0) = 

sign(tWo^o) X' where we restored the dimensional units. 
The Hall conductivity does not depend on the magnitude of 
the parameters t, Wq, and Aq, but depends on the sign of their 
product. 

When the pockets are present, the Fermi occupation func- 
tions in Eq. ( fTTI l exclude the pockets area from the integration 
over fc, which we represent symbolically as RBZ 

d^k B{k) 



, , f d^k 



4w^{k) ' 



(19) 



RBZ 



In this case, the Hall conductivity is no longer an integer 
At a finite temperature T, the Fermi occupation functions in 
Eq. ( fTTt further reduce the Hall conductivity. 

Now we turn to the calculation of the imaginary part of the 
ac Hall conductivity cr"y(ri). The imaginary part comes from 
the pole in Eq. ( fTTb at il = 2w{k). This pole represents verti- 
cal interband transitions induced by the electromagnetic wave. 
Taking the imaginary part of (fTTt . we find 



1 



J d^k B{k) 5[n-2w{k)]. (20) 



RBZ 



The integral over k in Eq. ( f20b is taken along the line in the 
Brillouin zone where the condition Vl — 2w{k) is satisfied. 
Eq. ^ is valid both for T = and T 7^ 0, as long as the 
condition :§> T is satisfied. 

Now we make a crude comparison of the theoretical result 
for the Kerr angle, given by Eqs. (|6]l and ( f20l l. with the exper- 
imental measurements [9]. We use the following representa- 
tive band structure parameters: t — 300 meV, t' — 90 meV, 
and /i = —260 meV. The refractive index is n « 1.69 lilST . 
and the interlayer distance is d = 1.17 nm. Assuming that 
the pseudogap originates from the ordered DDW state, we use 
the formula Wq{x) = 30(1 - x/xc) meV with Xc = 0.175, 
to model roughly how Wq varies with the hole concentration 
X fiod. The optical measurements were performed at the fre- 
quency nn = 0.8 eV. Taking Aq = r;Wo/2, with r; = 0.1%, 
we plot the value of the Kerr angle Ok vs. the hole doping x 
in Fig. [3] where we have multiplied the results by 2 to account 
for the contributions from two spin components. The magni- 
tude Ok roughly corresponds to the experimental values of the 
order of 100 nrad to 1 /irad for different dopings x ISJ. 
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FIG. 3: The calculated Kerr angle as a function of the hole concen- 
tration X. We use the linear dependence Wo{x) = Wo(l — x/xc) 
to model roughly the pseudogap variation with x, where Xc = 0.175 
and Wo = 30 meV. Although the admixture Ao = 0.001V|/o/2 of 
the dxy component is small, it leads to a non-zero Kerr angle of the 
same order of magnitude as in the experiments [9j. 



The intrinsic dc Hall effect predicted by Eq. ( fT9] l was not 
observed in the underdoped cuprates experimentally. We be- 
lieve this is because of the macroscopic domains with oppo- 
site chiralities present in a sample. Their contributions to the 
dc Hall effect would cancel out. On the other hand, the opti- 
cal measurements |9] are performed with a small laser beam 
within a single domain, so they give a non-zero effect. 

Within our theoretical model, the PKE is possible only 
when both the d^y and id^^-y'-i order parameters are present. 
Their transition temperatures are generally different, so the 
onset of the PKE corresponds to the lower of the two tran- 
sition temperatures. On the other hand, the higher transition 
temperature may be related to the "training" or "memory" ef- 
fect for the sign of the PKE [9]. The transition temperatures 
of the dxy and idr^2_y2 density waves should be taken as phe- 
nomenological parameters, because we cannot reliably calcu- 
late them from the microscopic first principles. 

In summary, we have developed a minimal theoretical 
model which can explain the macroscopic time-reversal- 



symmetry breaking leading to a non-zero polar Kerr effect 
||9|] in the underdoped cuprates. By detailed calculations, we 
show that admixing a very small component of the d^^j, order 
parameter into the postulated ordered DDW state picture of 
the pseudogap produces a PKE comparable in magnitude to 
that observed in Ref. |9]. Since the necessary dxy admixture 
is very small, only of the order of 0.1%, all other aspects of 
the YBCO phenomenology remain unaffected by our theory. 
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